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Abstract 


Since the molecular analysis of spike (S) glycoprotein gene of bovine coronavirus (BCoV) has been conducted and compared mainly 
among American and Canadian isolates and/or strains, it is unclear whether BCoV circulated in the other countries are distinctive in genetic 
characteristics. In the present study, we analyzed the S glycoprotein gene to characterize 10 winter dysentery (WD) coronavirus strains 
circulated in Korea during 2002—2003 and compared the nucleotide (nt) and deduced amino acid (aa) sequences with the other known BCoV. 
The phylogenetic analysis of the entire S glycoprotein gene revealed that the aa sequences of all Korean WD strains were more homologous to 
each other and were very closely related to respiratory bovine coronavirus (RBCV) strain OK and enteric bovine coronavirus (EBCV) strain 
LY-138, but were distinct from the other known BCoVs. Based on the phylogenetic analysis of the hypervariable region of the S1 subunit, all 
Korean WD strains clustered with the respiratory strain OK, BCQ3994 and the enteric strain LY-138, while the Canadian BCQ calf diarrhea 
and WD strains, and the American RBCV LSU, French EBCV F15 and avirulent VACC, L9, and Mebus strains clustered on a separate major 
branch. These data suggest that the WD strains circulated in Korea had a genetic property of both RBCV and EBCV and were significantly 
distinct from the ancestral enteric strain. 
© 2004 Elsevier B.V. All rights reserved. 
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Bovine coronavirus (BCoV), a member of the family Coro- 
naviridae, order Nidovirales (Cavanagh, 1995) is associated 
with severe diarrhea in newborn calves (CD), winter dysen- 
tery (WD) in adult cattle, and respiratory tract infections 
in calves and feedlot cattle (Cho et al., 2000; Clark, 1993; 
Lathrop et al., 2000; Saif et al., 1991; Storz et al., 
2000). BCoV possesses a single-stranded, enveloped, non- 
segmented RNA genome of positive polarity (De Vries et al., 
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1997). The virion contains five major structural proteins: the 
nucleocapsid (N) protein, the transmembrane (M) protein, 
the hemagglutinin/esterase (HE) protein, the spike (S) pro- 
tein, and the small membrane (E) protein (Lai and Cavanagh, 
1997). 

The variation in host range and tissue tropism of coron- 
aviruses is largely attributable to variations in the S glyco- 
protein (Gallagher and Buchmeier, 2001). The BCoV had 
point mutations or deletions especially in the S1 subunit of 
S glycoprotein (Chouljenko et al., 1998; Hasoksuz et al., 
2002; Gelinas et al., 2001; Rekik and Dea, 1994). Moreover, 
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molecular analysis of S gene of BCoV has been conducted 
and compared mainly between American and Canadian 
isolates and/or strains (Chouljenko et al., 1998; Hasoksuz 
et al., 2002; Gelinas et al., 2001; Rekik and Dea, 1994). In 
this study, we performed molecular analysis on the S gly- 
coprotein gene to characterize WD coronaviruses circulated 
in Korea during 2002-2003, compared the nucleotide, and 
predicted aa sequences with the other known BCoV. 

Ten WD. strains, KWD1-KWD10, were originally 
isolated in the G clone of human rectal tumor cells (HRT- 
18G) from fecal samples of adult dairy or beef cows, 
with WD in South Korea, and identified as a BCoV by 
immunoelectron microscopy, ELISA with BCoV-specific 
antisera and monoclonal antibodies, and RT-PCR specific 
for the part of BCoV N gene (Jeong et al., 2004). RNA 
was extracted from the supernatant fluids of HRT-18G 
cells infected with 10KWD strains, based on the acid 
quanidinium—phenol—chloroform RNA extraction method 
as described by Cho et al. (2001a). As a negative control, 
RNA was extracted from mock-infected HRT-18G cells. 
The oligonucleotide primers used in the RT-PCR were 
designed from the published sequence of the S gene of 
the Mebus strain (GenBank accession no. U00735). The 
sequence of primers (positions calculated from the start 
codon of the S gene) are as follows: Ist upstream primer, 
5'-ATGTTTTTGATACTTTTAATTTCC-3’ (1-24); 1st dow- 
nstream primer, 5’-ACACCAGTAGATGGTGCTAT-3’ (901— 
920); 2nd upstream primer, 5’-GGGTTACACCTCT- 
CACTTCT-3’ (782-801); 2nd downstream primer, 5’- 
GCAGGACAAGTGCCTATACC-3’ (1531-1550); 3rd ups- 
tream primer, 5‘-CTGTCCGTGTAAATTGGATG-3’ (1459 
1478); 3rd downstream primer, 5’-TGTAGAGTAATC- 
CACACGT-3’ (2268-2286); 4th upstream primer, 5’-TTC 
ACGACAGCTGCAACCTA-3’ (2151-2170); 4th down- 
stream primer, 5/-CCATGGTAACACCAATCC CA-3’ 
(3003-3022); 5th upstream primer, 5’-CCCTGT- 
ATTAGGTTGTTTAG-3’ (2691-2710); 5th downstre- 
am primer, 5’-ACCACTACCAGTGAACATCC-3’ (3587— 
3606); 6th upstream primer, 5’-GTGCAGAATGCT- 
CCATATGGT-3’ (3439-3459); 6th downstream primer, 
5'-TTAGTCGTCATGTGATGTTT-3’ (4073-4092). A 
one-step RT-PCR assay was performed as described by Cho 
et al. (2001a,b). The RT-PCR products were purified using 
a GenClean II kit (BIO 101, Inc., LaJolla, CA) according 
to the manufacturer’s instructions. The DNA sequencing 
was done using an automated DNA sequencer (ABI system 
3700, Applied Biosystem Inc., Foster City, CA). 

Using the DNA Basic module (DNAsis MAX, Alameda, 
CA), nucleotide sequences of our BCoV strains were 
first compared for the S sequence of BCoV as follows: 
the CD strains including Mebus (GenBank accession no. 
U00735), BCQ2070 (GenBank accession no. U06090), 
BCQ571 (GenBank accession no. U06093), BCQ20 (Gen- 
Bank accession no. U06092), BCQ9 (GenBank acces- 
sion no. U06091), and BCQ1523 (GenBank accession no. 
AAG40601); the WD strains including BCQ7373 (GenBank 


accession no. AAG40595) and BCQ2590 (GenBank acces- 
sion no. AKO1071); the RBCV strains including BCQ3994 
(GenBank accession no. AAK14398), LSU (GenBank ac- 
cession no. AAF25509), and OK (GenBank accession no. 
AAF25519), and EBCV strains including LY-138 (GenBank 
accession no. AAF25499) and F15 (GenBank accession no. 
D00731). The deduced aa sequences were then assembled 
and analyzed on the Amino Acid Basic module (DNAsis 
MAX, Alameda, CA). A sequence similarity search was per- 
formed for the BCoV S glycoprotein using the LALIGN 
Query program of the GENESTREAM network server at 
Institut de Génétique Humaine. Phylogenetic analyses were 
conducted using PhyloDraw program at the Graphics Appli- 
cation Laboratory, Pusan National University. 

All of the S genes of the 1OKWD strains contained an 
ORF of 4092 nucleotides, encoding a predicted protein of 
1363 aa residues and having a molecular weight of approxi- 
mately 150 kDa. Spike protein further divided it into $1 and 
S2 segments of approximately 86 and 65 kDa, respectively, 
at the cleavage site aa 768. Among all strains analyzed, a to- 
tal of 77 and 52 polymorphic nucleotides were identified in 
the S1 and S2 subunit of the spike glycoprotein of BCoV, re- 
spectively, compared with the Mebus strain (data not shown). 
These polymorphisms led to 46 and 24 aa changes at 41 and 
20 distinct sites, respectively. We have shown recently that the 
respiratory and enteric isolates collected concurrently from 
two different calves in USA were significantly distinct from 
the ancestral enteric strain, Mebus (Hasoksuz et al., 2002). 
By drawing comparisons of deduced aa between our strains 
and other known BCoV, the two most similar sequences were 
those of KWD10 strain and RBCV LSU strain (96.6%) in 
paired comparisons, and the most distant were those of the 
KWDS strain and the Mebus strain (93.4%). All the virulent- 
KWD strains tended to be distant from the ancestral enteric 
strain, Mebus, suggesting that allelic variation resulted in ge- 
netic drift over time (Hasoksuz et al., 2002). 

Based on the total number of aa substitutions, a phyloge- 
netic tree of all S gene sequences was analyzed with proto- 
type strain Mebus, RBCV strains LSU and OK, EBCV strains 
F15 and LY-138, and avirulent EBCV L9 (Fig. 1). The aa se- 
quences of 10 KWD strains were more homologous to each 
other and were very closely related to the virulent RBCV OK 
and EBCV LY-138 strains, among the other known BCoV 
studied in paired comparisons (Fig. 1). These data showed 
that the WD BCoVs circulated in Korea were genetically sim- 
ilar to both RBCV and EBCV. In addition, the phylogenetic 
comparison of Korean WD strains with other known BCoVs, 
on the sequences of the hypervariable region, revealed that 
all Korean WD strains were also clustered in groups with 
the virulent RBCV OK and BCQ3994 strains, and the viru- 
lent EBCV LY-139 strain (Fig. 2). The Canadian BCQ CD 
and WD strains, the American RBCV LSU strain, French 
EBCV F15 strain and avirulent VACC, L9 and Mebus strains 
were clustered on a separate major branch (Fig. 2). There- 
fore, these results indicated that the BCoV strains may be 
diverging from an enteric tropism to a dual (respiratory and 
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Mebus 


Fig. 1. The phylogenetic tree of the S gene of respiratory (RBCV), enteric (EBCV), calf diarrhea (CD), winter dysentery (WD), and avirulent strains or isolates 
was made using PhyloDraw at the Graphics Application Laboratory, Pusan National University. KWD1—KWD10: Korean WD strains; F15 and LY-138: EBCV 
strains; LSU and OK: RBCV strains; L9: avirulent vaccine strain; Mebus: prototype CD strain. 
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Fig. 2. The phylogenetic tree of the hypervariable region (aa 456-592; Rekik and Dea, 1994) of respiratory (RBCV), enteric (EBCV), calf diarrhea (CD), 
winter dysentery (WD), and avirulent strains. The CD strains including Mebus, BCQ2070, BCQ571, BCQ20, BCQ9, and BCQ1523; the WD strains including 
BCQ7373, BCQ2590, and KWD1-KWD10; the RBCV strains including BCQ3994, LSU, and OK; EBCV strains including LY-138 and F15. Phylogenetic 
analyses were conducted using PhyloDraw at the Graphics Application Laboratory, Pusan National University. 
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enteric) tropism over time via intermediates (Hasoksuz et al., 
2002). This speculation could be supported by our previous 
paper (Cho et al., 200 1a) in which gnotobiotic and colostrum- 
deprived calves were inoculated with respiratory isolates of 
BCoV shed viruses, both nasally and rectally. This subse- 
quently induced diarrhea. Moreover, feedlot cattle, which 
had shown respiratory disease signs at arrival time, devel- 
oped diarrhea 4 days after arrival, which peaked at 4-7 days 
post arrival (Cho et al., 2001b). Further study will be needed 
on the differences in pathogenicity, antigenicity, or any other 
physicochemical characteristics between the KWD strains. 
Inthe S glycoprotein, the hypervariable region spanning aa 
452-593 of the S1 subunit contains the S1B immunoreactive 
epitope which is the target for virus neutralizing monoclonal 
antibodies (Mabs) (Cavanagh, 1995). The hypervariable re- 
gion showed varied antigenicity when reacting with Mabs 
directed against the S glycoprotein. Therefore, Chouljenko 
et al. (1998) reported that four RBCV-specific aa substitutions 


at aa 510, 531, 543, and 578 can be predictive and Mabs- 
specific strains for this region may also be distinguished 
among respiratory, enteric, and vaccine BCoV strains. In 
the present study, three aa sites including aa 510, 543, and 
578 of RBCV OK and LSU strains were identical in all Ko- 
rean WD strains (Fig. 3). In our previous paper (Hasoksuz 
et al., 2002), enteric and respiratory BCoV pairs collected 
concurrently from two feedlot calves in USA also had iden- 
tical aa substitution sites at the same sites. Furthermore, aa 
531 (N — G) substitution of RBCV OK and LSU strains 
compared with the Mebus strain was also detected in the CD 
strain BCQ 1523 (data not shown). Since these changes were 
observed in WD, CD, and RBCV strains, our data does not 
support the speculation of Chouljenko et al. (1998). In addi- 
tion, aa 13 and 744 of the S1 subunit have been reported for 
EBCV-specific sites, which were observed in EBCV strains 
F15 and LY (Chouljenko et al., 1998). However, aa 13 and 
744 of the S1 subunit of all Korean WD strains were identical 
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Fig. 3. Comparison of the predicted amino acid sequences of the BCoV S glycoprotein specified by different strains. Light-gray boxes contain aa sites known 
for respiratory bovine coronavirus-specific (Chouljenko et al., 1998), dark boxes contain virulent-specific (Chouljenko et al., 1998), and clear boxes contain 
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with those of LSU and OK strains of RBCV, BCQ7373, and 
BCQ2590 of WD isolates, Mebus, and L9 avirulent strains 
(Fig. 3). In particular, aa 744, reporting its sequence in many 
isolates and/or strains, are the same as other known BCoV 
including RBCV, CD, WD, and EBCYV. These results imply 
that these aa sites are EBCV FI15 and LY strain-specific, or 
that there are no EBCV-specific sites in the S gene. 

The N-terminal region of the S 1 subunit of 10 KWD strains 
(aa 1-330; Laude et al., 1995), being shown to function as a 
receptor-binding domain in mouse hepatitis virus, had a total 
of 15 aa changes compared with the Mebus strain (data not 
shown). Amino acid substitutions in this region could alter 
the tropism of the coronavirus. For example, the first 200 aa 
of the S1 subunit of the PRCV were associated with respira- 
tory tropism and reduced enteropathogenicity (Wesley et al., 
1991). Therefore, these 15 aa substitutions in the Korean WD 
strains compared with that of the Mebus strain may induce 
either an alteration of the receptor-binding ability during viral 
invasion or pathogenicity to cattle. 

The S protein of BCoV has a cleavage site located at aa 
763-768 (KRRSRR) and is divided into the peripheral S1 
glycoprotein from the transmembrane domain of the S2 gly- 
coprotein (Abraham et al., 1990). The sequence KRRSRR at 
the predicted proteolytic cleavage site was identified in all 
KWD strains. Interestingly, all KWD strains had an Ala to 
Ser change at aa 769, immediately upstream of the cleav- 
age site. Since this change was observed exclusively in the 
respiratory isolates LSU and OK, it was speculated that the 
change in aa 769 was a potential marker of respiratory tropism 
(Chouljenko et al., 1998). Although changes in aa surround- 
ing the cleavage site of the HA molecule of the influenza A 
virus and Sendai virus have been demonstrated to modulate 
its cleavability (Rott et al., 1984; Hsu et al., 1987), Rekik and 
Dea (1994) reported that this change (769; A — S) did not 
appear to modulate cleavability, the rate of viral replication 
or the type of cytopathic changes induced in HRT-18 cells for 
their enteric BCoV isolates. In our previous paper (Hasoksuz 
et al., 2002), this change (769; A — S) was observed in both 
respiratory and enteric strains. Taken all together, the aa sub- 
stitution (A — S) at aa 769 does not appear to be a potential 
marker of respiratory tropism. 

Of the two cleavage products of S, the $2 subunit is highly 
conserved among coronaviruses (Gallagher and Buchmeier, 
2001). In the present study, the aa sequences of the S2 sub- 
unit of Korean WD strains were comparatively well con- 
served compared with the $1 subunit of the Mebus strain. 
The notable aa changes were observed only in the first hy- 
drophobic region of all KWD strains. In this region, there 
were four aa substitutions (aa 959, 965, 984, and 988), re- 
sulting in a marked increase of hydrophilicity in comparison 
with that of the Mebus strain (data not shown). Interestingly, 
the other known virulent BCoV, including RBCV LSU and 
OK strains, EBCV LY and F15 strains revealed an increase 
in hydrophilicity in the first hydrophobic region compared 
with avirulent strains Vacc, L9, and Mebus (data not shown). 
The membrane fusion activity is most likely conferred by the 


internal hydrophobic sequences within the S2 subunit (Luo 
and Weiss, 1998; Yoo et al., 1991). Therefore, aa substitu- 
tions inducing an increase in hydrophilicity observed in our 
strains and other virulent BCoV may alter the fusion activity 
or pathogenicity of the viruses. 

Putative virulent-specific sites of BCoV were reported 
to exist at 10 aa sites (six within the S1 subunit and 
four within the S2 subunit) in the S gene through com- 
parison with avirulent and virulent strains (Zhang et al., 
1991). In the present study, the S glycoprotein of all KWD 
strains contained seven aa substitutions (aa 33, 40, 248, 470, 
965, 1241, and 1341) which were common for all virulent 
strains (Fig. 3). This is in agreement with Chouljenko et al. 
(1998), who identified seven virulent-specific aa changes in 
the S gene by comparison with the RBCV, EBCYV, and vac- 
cine strains. Since these virulent-specific sites were observed 
in all virulent BCoV, but not in avirulent BCoV, these aa 
changes may affect the structure and function of the S glyco- 
protein and alter the pathogenetic potential of these viruses 
(Chouljenko et al., 1998; Zhang et al., 1991). 
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